Gamma-rays from the compact colliding wind 
region in Cyg OB2 #5 

A. T. Araudo, G. N. Ortiz-Leon and L. F. Rodriguez 

Centra de Radioastronomia y Astrofisica, Universidad National Autonoma de Mexico, A.P. 3-72 
(Xangari), 58089 Morelia, Michoacdn, Mexico 

Abstract. In this contribution we model the non-thermal emission (from radio to y-rays) pro- 
duced in the compact (and recently detected) colliding wind region in the multiple stellar system 
Cyg OB2 #5. We focus our study on the detectability of the produced y-rays. 
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INTRODUCTION 

The collision between stellar winds of massive stars is a promising scenario to produce 
non-thermal emission. In early type binaries, a strong shock is created between the 
stars, where the wind ram pressures are equated. In those shocks particles can be 
accelerated up to relativistic energies via the Fermi-I acceleration mechanism, and 
produce non-thermal emission. In particular, y-rays can be produced by the interaction 
of the accelerated particles with the photon and matter fields provided by the stellar 
radiation and wind, respectively. 

Non-thermal radio emission from colliding wind regions (CWRs) has been firmly 
detected from many systems and they are putative candidates to be y-ray sources since 
many years ago. In the past, only marginal correlations between CWRs and unidentified 
y-ray sources (of e.g. EGRET catalog) has been found. However, the recent detection of 
Eta-Car by Fermi [1] and Agile [2] shows that the y-ray emission of CWRs is not only a 
theoretical prediction. In this work we study other promising y-ray source to be detected 
by Fermi: the peculiar system Cyg OB2 #5. 



The source Cyg OB2 #5 

Cyg OB2 #5 is a radio-bright early-type multiple system located into the big stellar 
association Cyg OB2, at a distance of ~ 1.4 kpc. This peculiar source is a (possible) 
quadruple stellar system that shows radio structures over different scales, from about 10 
milliarcseconds (mas) to about 30 arcseconds [3]. At the 10 mas scale, there is a non- 
thermal, arc-like structure that traces the CWR between the wind of an eclipsing contact 
binary and that of an unseen nearby companion. About \" to the NE of the contact binary, 
there is another non-thermal component that results from the interaction of the winds of 




FIGURE 1. Left: Compact and extended components of the system Cyg OB2 #5. The sketch is not 
drawn to scale. Right: Escape, radiative cooling and acceleration times for the cases discussed in the text. 



the contact binary and that of a known B-type star. Finally, there is an extended (~ 30") 
synchrotron structure that could also be associated with Cyg OB2 #5. In Figure 1 (left) 
we show a sketch of the three non-thermal structures that compose this complex source. 
In this study we model the y-ray emission produced in the compact CWR (produced by 
the contact binary and an unseen companion). 



PARTICLE ACCELERATION AND LOSSES 

Both stars of the compact CWR are separated a distance D = x v + x s , and the CWR is 
located at a distance x p and x v = a//3 x s from the primary and the secondary components 
of the system, respectively. Being /3 = MpVoo.p/MsVoo.s, where M is the mass-loss rate, 
and Voo is the terminal wind velocity, the stagnation point is located closer to the less 
powerful star (the unseen companion in our case). Characterizing the wind of the contact 
binary by v^.p ~ 1500 km s _1 and M p ~2x 10~ 5 M yr _1 , and considering that the 
companion is a BO. 5 early-type star [4] x s results ~ 1.5 AU, where we have estimated 
x p = 25 AU from radio data. At the location of the stagnation point, a discontinuity 
surface and a double bow-shock structure is created. In long period binary systems, the 
formed CWR is thick and it is not affected by disruptive instabilities. Thus, CWRs are 
good places for particle acceleration, since at least one of the shocks is adiabatic and the 
shocked region will be sufficiently large to allow charged particles to reach relativistic 
energies. In this work we assume Fermi-I as the operating mechanism to accelerate 
particles. Being the wind of the contact binary the most powerful, we will consider only 
the acceleration of particles in the shock produced by this wind. 

The detection of synchrotron emission is an evidence that there is a population of 
relativistic electrons (e) in the CWR. However, protons (p) can be accelerated as well. 
Particles accelerated by the Fermi-I mechanism follow a power-law energy distribution 
Ne :P = K e>p E~p , where we have assumed that both electrons and protons are injected 
with the same spectral index 2.1. Considering that the same number of electrons and 
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FIGURE 2. Spectral energy distributions for the cases discussed in the text. The Fermi sensitivity is 
also plotted in order to show that the source is detectable by this instrument. 



protons are accelerated at energies < m e c 2 [5], and neglecting ionization/coulomb losses, 
we relate K p = {m p /m e )^ 2A -^l 2 K e . 

From the measured synchrotron flux S v (v = 8.4GHz) = 3 mJy [4] we estimate the 
magnetic field B in the CWR by equating magnetic and non-thermal energy densities: 
Ub = U e + aUp, where U e and U p are the energy densities of relativistic electrons and 
protons, respectively, and a takes two values: (no protons are accelerated) and 1 
(protons are accelerated as well as electrons). We obtain B = 0.07 G in the former 
case and 0.29 G in the latter. Once we know B, we can calculate the acceleration and 
synchrotron cooling time. To calculate the losses by Inverse Compton (IC) scattering, 
the largest density of target photons (£/ pn = 4.45 erg cm~ 3 at the location of the CWR) 
is provided by the unseen companion with luminosity L s ~ 6. 1 x 10 4 L [4]. In addition 
to synchrotron and IC losses, accelerated electrons in the field of ions in the wind can 
produce relativistic Bremsstrahlung radiation at y-ray energies. However, the density 
of these ions in the CWR is not enough (n w = 1.3 x 10 7 cm~ 3 ) for the relativistic 
Bremsstrahlung losses to be more efficient than IC cooling. Thus, the maximum energy 
of accelerated electrons is constrained by IC losses, given E e , max ~ 40 and 93 GeV in 
the cases with a = and 1, respectively, as is shown in Fig. 1 (right). 

Relativistic protons interacting with ions in the stellar winds produce neutral pions 
and then y-ray s. However, relativistic hadrons escape from the CWR on a time shorter 
than the pp cooling time. At energies E p < 5 TeV, advection losses are dominant 
fedv ~ 4x p /voo.p ~ 10 7 s), and at larger energies the diffusion escape become more 
relevant. This produces that N p shows a break at energy = 5 TeV. The maximum 
energy of protons is constrained by diffusion losses given E p max ^ 

50 TeV. 

GAMMA-RAY EMISSION AND DISCUSSION 

We computed the spectral energy distribution (SED) of synchrotron, IC, relativistic 
Bremsstrahlung, and pp, in the CWR considering the distributions N e and N p derived 



previously by adjusting the synchrotron spectrum with radio data and below the equipar- 
tition condition. We obtain that the most important radiation mechanism to produce y- 
rays in the compact CWR of Cyg OB2 #5 is the IC scattering due to the dense photon 
stellar field provided by the unseen companion when only electrons are considered (case 
with a = 0). In this case synchrotron, IC and relativistic Bremsstrahlung bolometric lu- 
minosities result L S y n ~ 1.5 x 10 30 , Lie ~ 2.6 x 10 34 , and L B rem ~ 5.4 x 10 31 erg s 
respectively, as is shown in Fig. 2 (left). In the case with a = 1, the non-thermal en- 
ergy is distributed between electrons and protons resulting that U p > U e as a conse- 
quence of m p > m e [5]. The achieved bolometric luminosities are L syn ~ 1.7 x 10 30 , 
Lie ~ 1-7 x 10 33 , and L^ Km ~ 2.9 x 10 30 erg s _1 . In this case proton-proton collisions 
are the most important cooling channel in the y domain, reaching a bolometric luminos- 
ity L pp ~ 7.8 x 10 32 erg s _1 , as is shown in Fig. 2 (right). 

The y-ray spectrum is attenuated at energies greater than ~ 0.3 TeV, where the optical 
depth due to photon-photon electron-positron pair production (y + y — y e + + e~) in the 
intense stellar radiation field of the BO. 5 companion star reaches non-negligible value. 
In the case with a = 0, photon-photon absorption is not relevant because E ph ~ 0.3 TeV 
is larger than the maximum energy of photons produced there, but in the case with a = 1 
the attenuation is significant. 

From Fig. 2 it can be seen that the y-ray emission produced in the compact CWR of 
Cyg 0B2 #5 (by the simple model described in the present contribution) can be detected 
with the Fermi satellite. In particular, when only the relativistic electrons (that produce 
the detected synchrotron emission) are considered, IC scattering of photons from the 
unseen companion reach specific luminosity values in the GeV domain detectable with 
Fermi in a deep exposure. When protons are also accelerated up to relativistic energies, 
y-rays from pp collisions are also detectable by Fermi. These results shown that even 
in the case in which the Fermi detection from the Cyg OB2 region has been associated 
with the pulsar PSR J2032+4127, under certain conditions a significant amount of y- 
rays can be produced in Cyg OB2 #5. The detection of y rays from CWRs can help us 
to understand the properties of stellar winds, as well as the development of high-energy 
processes in stellar environments. 

ACKNOWLEDGMENTS 

This work is supported by CONACyT, Mexico, and PAPIIT, UNAM. 

REFERENCES 

1. Abdo et al. [Fermi collaboration] 2010, ApJS 188, 405 

2. Tavani, M, Sabatini, S., Pian, E., et al. 2009, ApJL, 698, L142 

3. Ortiz-Leon G.N., Loinard L., Rodriguez L.F., Mioduszewski, A.J. & Dzib A.S. 201 1, ApJ 730, 30 

4. Dzib, A.S., et al. 2012 (submitted) 

5. Bell, A.R. 1978, MNRAS 182, 443 



